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Employing the trace distance as a measure for the distinguishability of quantum states, we study 
the influence of initial correlations on the dynamics of open systems. We concentrate on the Jaynes- 
Cummings model for which the knowledge of the exact joint dynamics of system and reservoir allows 
the treatment of initial states with arbitrary correlations. As a measure for the correlations in the 
initial state we consider the trace distance between the system-environment state and the product 
of its marginal states. In particular, we examine the correlations contained in the thermal equilib- 
rium state for the total system, analyze their dependence on the temperature and on the coupling 
strength, and demonstrate their connection to the entanglement properties of the eigenstates of the 
Hamiltonian. A detailed study of the time dependence of the distinguishability of the open system 
states evolving from the thermal equilibrium state and its corresponding uncorrelated product state 
shows that the open system dynamically uncovers typical features of the initial correlations. 

PACS numbers: 03.65.Yz,03.65.Ta,42.50.Lc 



I. INTRODUCTION 

Quantum systems are typically subjected to the inter- 
action with an environment which influences their dy- 
namics in a non negligible way. A realistic descrip- 
tion, taking this external influence into account, is cru- 
cial for the theoretical description of open quantum sys- 
tems, which play an important role in many areas of 
physics pp. The available theoretical tools allow a full 
characterization for a Markovian dynamics, which can 
be described by means of completely positive quantum 
dynamical semigroups However, the assumptions 

justifying the Markovian description of the open system 
dynamics are often too restrictive, and a more general 
analysis is required. A wealth of different approaches to 
deal with non-Markovian dynamics have been introduced 
but they typically rely on the hypothesis that at 
the initial time the open system and the environment are 
statistically independent. This assumption is well justi- 
fied in the case of weak interaction, but in general one 
cannot neglect the initial correlations between the open 
system and the environment ffifl426| . 

Considering an open system S which is coupled to an 
environment E and assuming that the composite system 
evolves according to a unitary time evolution operator 
U(t) from a total initial state pse(0), we can write the 
state of S at time t as follows, 

Ps (t) = Tr E [U(t)p SE (0)UHt)}. (1) 
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This equation defines a linear, completely positive and 
trace preserving map A t from the state space of the total 
system S + E to the state space of the open system S : 

Pse{0)^ p s {t) = A tPSE {0). (2) 

If the open system S and its environment E are initially 
in an uncorrelated tensor product state 

Pse(0) = Ps (0)®pe (3) 

with a fixed environmental state p E , Eq. ([!]) also defines 
a linear map <I> t from the state space of S into itself, 

Ps(0) ^ Ps(t) = $tPs{0) = Tr B [U(t)p s (p) ® p E UHt)} . 

(4) 

It can be shown that this quantum dynamical map 4> t 
is again completely positive and trace preserving. Un- 
der the additional assumption that the family of dynam- 
ical maps > 0} constitutes a semigroup, one de- 
rives the general mathematical structure of its generator, 
which leads to the widely-used quantum Markovian mas- 
ter equations for the open system state ps(t) in Lindblad 
form. 

The above construction of the quantum dynamical 
map $t presupposes that one restricts the class of ini- 
tial conditions to states of the form of Eq. ([3|, where pe 
is a fixed, given environmental state. Thus, a large class 
of initial conditions is excluded when considering dynam- 
ical maps acting on the reduced state space, in particu- 
lar those initial conditions that describe correlations and 
entanglement between system and environment. On the 
other hand, it is a well known fact that correlations in the 
initial state can have strong influences on the open sys- 
tem dynamics, both in thermal equilibrium and in non- 
equilibrium systems. The question is thus, how do ini- 
tial correlations affect the reduced system dynamics, and 
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what are appropriate observable measures that quantify 
such effects? Here, we discuss these questions in detail 
with the help of the example of the Jaynes-Cummings 
model, the model of a two-state system coupled to a 
bosonic field mode, employing the analytical representa- 
tion of the reduced dynamics of this model for arbitrary 
initial states [2"7j . 

It should be mentioned that under certain additional 
assumptions Eq. (JlJ can indeed be used to construct 
maps on the reduced state space to represent the dy- 
namics in the case of initial correlations |28H33] . How- 
ever, this construction demands that the initial correla- 
tions between the system and its environment are fixed. 
It turns out that the dynamical maps arising in this way 
may be not completely positive, and not even positive. 
This requires the determination of a certain compatibil- 
ity domain in the physical state space, which is a very 
complicated mathematical task. In the present paper we 
shall follow an entirely different strategy, to analyze the 
role of initial correlations. Namely, in order to quantify 
the effect of initial system-environment correlations in 
the subsequent time evolution of the open system, we will 
investigate the trace distance D{p 1 s {t) ) p 2 s {t)) between a 
pair of states pg(t) and pg(t) of S, which evolve from a 
given pair of initial states p\ E (0) and p 2 SE (0) of the total 
system. This approach has also been used in [TTJ [TB] to 
construct a measure for the non-Markovianity of quan- 
tum processes, and in [34 to develop a witness which 
allows the detection of initial correlations through only 
local measurements on the open system. An application 
to a specific system has been recently considered in |35j . 

In the present paper we will address, in particular, 
the situation in which pg E (0) represents a thermal equi- 
librium (Gibbs) state corresponding to the full Hamil- 
tonian of the model. Since this correlated state is 
invariant under the time evolution, its reduced states 
pUO) = TW^(O) and ^(0) = Tr s p^(0) remain 
of course time independent. However, the initial state 
Pse(®) = Psi®) ® Pe(®) gi yen by the product of the 
marginals does evolve in time, and we will investigate 
the dynamics of the trace distance D(pg(t), p%(t)) be- 
tween the open system states corresponding to the ini- 
tial states Pse(Q) an( i Pse(O)- For this case the trace 
distance is bounded from above by the trace distance 
D(p 1 SE (0), pg(0) ® Pb(0)) which provides a measure for 
the amount of correlations in the initial Gibbs state [51] . 
Analyzing in detail the dependence on the temperature 
and on the system-environment coupling strength, we 
demonstrate that at small temperatures characteristic 
properties of these correlations are related to the eigen- 
value spectrum and, in particular, to the quantum corre- 
lations and the entanglement structure of the eigenstates 
of H . We will discuss further the signatures of these 
properties in the subsequent dynamics of the open sys- 
tem states. It will be shown that, in fact, the open system 
dynamically uncovers typical features of the correlations 
in the initial states. 

The paper is organized as follows. In Sec. [il] we intro- 



duce the trace distance and show its relevance as a mea- 
sure of the distinguishability of two quantum states and 
of the correlations contained in a given bipartite state. 
We further consider the exact reduced dynamics of the 
Jaynes-Cummings model, and study as an example the 
time behavior of the distinguishability of distinct initial 
states. In Sec. |III| we provide a detailed study of the 
correlations contained in the Gibbs state associated to 
the Jaynes-Cummings Hamiltonian, as measured by the 
trace distance between the state and the tensor product 
of its marginals. We then investigate the time depen- 
dence of the distinguishability of the corresponding time 
evolved states. 



II. TRACE DISTANCE AND INITIAL 
SYSTEM-ENVIRONMENT CORRELATIONS 

A. General theory 

1. Properties and physical interpretation of the trace 
distance 

The trace distance of two trace class operators A and 
B is defined as \ times the trace norm of A — B, 

D(A,B) = \\\A-B\\ U (5) 

where the trace norm of an operator X is defined by 

||X||i = Tr|A| = TWX^X. (6) 

If X is trace class and self-adjoint with eigenvalues Xi, 
this formula reduces to the sum of the absolute eigenval- 
ues (counting multiplicity), 

llAUi =X>I- (7) 

i 

The trace distance of two quantum states, represented 
by positive operators p 1 and p 2 with unit trace, is thus 
given by 

D(p\p 2 ) = ±\\pi - p 2 \\ 1 = ±Tr\pi - p% (8) 

The trace distance is a metric on the space of physical 
states with several nice properties which make it a useful 
measure for the distance between two quantum states. 
We list some of them: 

1.: The trace distance for any pair of states satisfies the 
inequality 

0<D(p\p 2 )<l, (9) 

where D(p ,p ) = if and only if p 1 = p 2 , and 
D(p 1 , p 2 ) — 1 if and only if p 1 and p 2 have orthog- 
onal supports. 
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2.: Being a metric, the trace distance satisfies the trian- 
gular inequality, 



D{p\p*)<D(p\p A )+D[p\p 2 ) 



(10) 



3.: All trace preserving positive maps A are contractions 
of the trace distance 361, 



D(A.p x ,Apr) < D{p L ,p 2 ), 



(11) 



where the equality sign holds if A is a unitary trans- 
formation. 

4. : The trace distance is subadditive with respect to the 

tensor product, 

Dip 1 ®a\p 2 ®a 2 ) <D(p\p 2 )+D{a\a 2 ). (12) 
In particular, one has 

D(p 1 ®a,p 2 ®a) = D{p 1 ,p 2 ). (13) 

5. : The trace distance can be represented as a maximum 

taken over all projection operators II, 

D(p\p 2 )=m^Tr{U(p l -p 2 )}. (14) 

The physical interpretation [57] of the trace distance 
is based on the relation (14 1. Suppose Alice prepares 
a system in one of two quantum state p 1 and p 2 with 
probability of 1/2 each. She gives the system to Bob, 
who performs a measurement in order to distinguish the 
two states. Employing Eq. (14) one can show that the 



maximal success probability for Bob to identify correctly 
the state is given by [l + D(p 1 , p 2 )] /2. This means that 
the trace distance represents the maximal bias in favor 
of the correct state identification which Bob can achieve 
through an optimal strategy. Hence, the trace distance 
D{pi,p2) can be interpreted as a measure for the distin- 
guishability of the states p 1 and p 2 . 



2. Dynamics of the trace distance 

We consider any two total system initial states p$ E (0) 
and Ps E {0), and the corresponding open system states 
Ps(t) and Ps(t) at time t. According to Eqs. ([!]) and ^ 
the latter are given by p\{t) = Ktp l SE {Q) and p|(t) = 
A tj o| B (0). Since A t is completely positive and trace pre- 
serving, we obtain from Eq. (11) a bound for the trace 



distance between the reduced system states, 

D(J s {t), p% (*)) < D(pl E (Q), p 2 se (0)). (15) 

If the initial states are uncorrelated with the same envi- 
ronmental state pe, that is P5 E (0) = pg(0) <g> Pe(0) and 
Pse(®) = ® Ae(0), this inequality reduces with the 



help of ( 13 1 to the contraction property for the dynamical 
map Q, 



This means that for initially uncorrelated total system 
states and identical environmental states the trace dis- 
tance between the reduced system states at time t can 
never be larger than its initial value. 
The inequality ( 15 ) may be written as 



£(4(t),^(t))-D(^(0),p|(0)) 

< D(p^(0), p 2 SE (0)) - £(Ps(0), Ps(0)) 

= /(pk(0),p| E (0)). 



(17) 



According to this inequality the change of the trace dis- 
tance of the open system states is bounded from above by 
the quantity I(ps E (0), p| B (0)) > 0. This quantity rep- 
resents the distinguishability of the total system initial 
states minus the distinguishability of the corresponding 
reduced system initial states. Thus, I(pseiPse) can t> e 
interpreted as the relative information of the total ini- 
tial states which is initially outside the open system, i.e., 
which is inaccessible for local measurement performed on 
the open system [54] . 

For I(pg E (0), Pg E (0)) > the trace distance of the 
open system states can increase over its initial value. 
This increase can be interpreted by saying that infor- 
mation which is initially outside the open system flows 
back to the system and becomes accessible through lo- 
cal measurements. Note that, as will be illustrated by 
means of several examples below, the bound for the dy- 
namics of the trace distance given by Eq. (17) is tight, 
i.e., it can be reached for certain total initial states. If the 



bound of inequality ( 17 1 is actually reached at some time 



t, the initial distinguishability of the total system states is 
equal to the distinguishability of the open system states 
at time t. This means that the relative information on 
the total initial states has been dynamically transferred 
completely to the open system [34] . 



Using the sub-additivity of the trace distance ( 12 ) and 
the triangular inequality (10) one deduces from (17) the 
following inequality [34] . 

^(PsW,P5W)-^(Ps(0),p|(0)) 
<D(p l SE (0),pl(0)®p E (0)) 

+D(p 2 SE (0), P 2 (Q) ® p%(0)) + D{p E {0),p 2 E {0)). (18) 

For any state pse the quantity D(pse, Ps®Pe) describes 
how well P se can be distinguished from the fully uncorre- 
lated product state P s® Pe of its marginal states ps and 
Pe- Thus, D(pse, Ps ® Pe) can be interpreted as a mea- 
sure for the total amount of correlations in the state P se- 
Therefore, the inequality ( 18 1 shows that an increase of 



the trace distance of the open system states over its ini- 
tial value implies that there must be correlations in the 
initial states Pse(®) or or t na t, the environmen- 

tal states are different. An important special case, which 
will be considered in detail in the present paper, occurs 
if p 2 SE (0) is given by the product state obtained from the 
marginals of p l SE {Q), i.e., p 2 SE (0) = ,4(0) <g> p E (0). The 
inequality (|17|) then reduces to the simple form 



D(pUt),p s (t)) < D(pl(0),p 2 s (0)). 



(16) 



D{ Ps (t), P 2 s {t)) < DG^(0), pUO) ® Pe(0)), (19) 
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according to which the increase of the trace distance is 
bounded by the amount of correlations in the total initial 
state [Ml. 



B. Example: The Jaynes-Cummings model 

1. The physical model 

We consider a two-state system coupled to a single 
mode of the radiation field with total Hamiltonian 

H = Hs + H E + Hj 

= uj <t + o-- + w6 t & + g(cr + <g)fe + er_ ® 6 f ) , (20) 

where cr+ = |1)(0| and er_ = |0)(1| are the raising and 
lowering operators of the two-state system, b' and b are 
the creation and annihilation operators of the field mode, 
and the coupling term is in the Jaynes-Cummings form. 
This model describes, e.g., the interaction between a two- 
level atom and a mode of the radiation field in the electric 
dipole and rotating wave approximation. In the interac- 
tion picture the Hamiltonian takes the form 



Hi{t) = g (<7+ ® be iAt + cr_ <g ftt e -<A*) 



(21) 



where A = loq — uj denotes the detuning between the 
system's transition frequency luq and the frequency ui of 
the field mode. The exact time-evolution operator for 
the total system in the interaction picture can then be 
written as (see, e.g., Ref. [55]): 



U(t) = 



c(n+l,t) d(n + l,t)b 



(22) 



where we have introduced the following functions of the 
number operator h = wb, 



c(n, t) 
d{h, t) 
with 



_ AAt/2 



cos I O(ri) ; 



i — , „ ; sin [ fi(n) — 



fi(n) V 2 



(23) 



Q(h) = VA 2 + Ag 2 n. (24) 
With the help of the unitary time-evolution operator 



given by Eq. ( 22 1 we can easily determine the exact ex- 



pression for the reduced density matrix of the two-level 
system at time t, 



Ps{t) 



Pu(t) pw(t) 
P*o(*) Poo(t) 



(25) 



corresponding to an arbitrary initial state pse(0) of the 
total system. First, we expand ps_e(0) with respect to 
the basis vectors |a) ® |n) = \a, n), where a = 1, labels 



the states of the two-state system, and n = 0, 1, 2, . . . the 
number states of the field mode, 

Pse(0)= J2 Pa™(0)|a,m)(/3,n|. (26) 

Substituting this expression into Eq. ([T]) with U (t) given 
by Eq. (p2|), one obtains 



Pn(t) = E (0)|c»+i(t)| 2 + 2vVTT 
n 

x Re {p^ n+ \0)d* n+1 (t)c n+1 (t)} + np^{Q)\dn(f)f 
Pw{t) = [-vVTT i o" 1 +1 '"(0)c rl+2 (i)d„ + i(t) 

n 

-V^^V^+iPoi 2 ' n (o)d n+2 (t)d n+1 {t) 



(27) 

where c n (t) and d„(<) denote the eigenvalues of c(n,t) 
and d(n,t) corresponding to the eigenstate |n), respec- 
tively. 



We note that Eq. (271 does in general not lead to a 



dynamical map for the state changes of the reduced two- 
state system since it is not possible to write the right- 
hand side of this equation as a function of the matrix 
elements of the reduced initial state ps(0) which are given 
by 



A 



a/j(o) = 5>S2(o). 



(28) 



However, if the total initial state is of tensor product 
form, pse(0) = Ps(0) <S> Pe(0) and, therefore, 



Pa™(0)=PaM0)p" m (0), 



(29) 



it is indeed possible to construct the dynamical map; if 
moreover [p^(0),n] = 0, one finds the map already de- 
rived in Ref. [2T] . 



2. Dynamics of the trace distance for pure or product total 
initial states 



We illustrate the dynamics of the trace distance and 
the inequality ( 17 ) by means of two simple examples, con- 



sidering the situation in which the total initial state is a 
product state or a pure state. The case of a mixed, corre- 
lated initial state will be considered in detail in Sec. IIIIl 



The quantity on the right-hand side of Eq. (171, rep- 



resenting the information which is initially outside the 
reduced system, can be larger than zero basically for two 
reasons: First, because one has different environmental 
initial states p E (0) and p E (0) and, second, because of 
the presence of correlations in the initial states pg E (0) 
or Pg E (0) (see inequality (|l8|). To illustrate the first 
case we study the trace distance between the two reduced 
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states pg(t) and p 2 s (t) evolving from two product initial 
states with the same reduced system state, namely from 
PseW = Ps(0) ® ^(0) and p 2 SE (0) = Ps (0) ® p|(0), 
where 

Ps (0) = |a 1 | 2 |0)(0| + |/3 1 | 2 |1)(1| (30) 

and the two environmental states are taken to be 

p E (0) = k| 2 |n)(n| + |A| 2 |« - l)<n- 1|, i = 1,2, 

(31) 

with the normalization condition \ai\ 2 + \f3i\ 2 — 1. Nu- 
merical simulation results for this case are shown in 
Fig. (flja). We see from the figure that the bound of 
Eq. |l7L which is given by ||ai| 2 — |a2| 2 |, is indeed 
reached here. For a study of the second case we consider 
an initially correlated pure state of the form 

Pse(0) = MM. (32) 

with \tp) =a\Q,n)+p\l,n-l), \a\ 2 + \/3\ 2 = 1, together 
with an initial product state of the form 



p| E (0)=p|(0)®p|(0) 



(33) 



|/3| 2 |0)(0| 



1)1 



and 

^2 



P 2 e(0) 



with p|(0) 

|a| 2 |n)(n| + \/3\ 2 \n - l)(n - 1|. Note that p 2 SE (0) is not 
equal to the product of the marginals of p 1 SE (0) . As can 
be s een from Fig. ((TJb) also for this case the bound of 



Eq. (17), which is given by 1(1 + \a\ A + |/3| 4 ), is repeat- 



edly reached in the course of time. As expected, in both 
cases the trace distance of the states exceeds its initial 
value, corresponding to the fact that the reduced system 
dynamically retrieves the information initially not acces- 
sible to it, related to the different initial environmental 
states or to the initial system-environment correlations. 
Note that the trace distance starts increasing already at 
the initial time, indicating that the information is flow- 
ing to the reduced system from the very beginning of the 
dynamics. Moreover, it keeps oscillating also for large 
values of t, so that the distinguishability growth between 
reduced states can be detected, e.g. by quantum state 
tomography, also making observations after a long inter- 
action time. 

In both situations considered and visualized in Fig. ([!]) 
the maximum value of the trace distance as a function 
of time is equal to the upper bound given by Eq. ( fl7| , 
indicating that the information initially inaccessible to 
the reduced system has been transferred completely to 
it during the subsequent dynamics. This is of course 
not always the case and it is an important problem to 
characterize explicitly those initial states for which such 
a behavior indeed occurs. Let us consider the special 
case given by Eq. (19 1, in which the two total initial 



states are a correlated state and the tensor product of its 
marginals, taking p l SE {0) to be a pure entangled state, 
i.e., Pg E (0) = \ifj)(ip\ with = a\0,n) +/3|l,m). For 
this case Eq. (27) leads to 



D(p 1 s (t) 1 p 2 s (t)) = 
-\c n+1 {t)\ 2 ) + 25 n 



|a/3| 2 (| Cm+1 (i)| 2 ~| C „(i)| 2 4 
yn -i^/nRe{a*/3d* n (t)c n (t)}\ 



m 2 

(34) 



D( P ;(i). P |(t)) 

0.10 
0.08 
0.06 
0.04 
0.02 



D(pi(l).p|(t)) 





FIG. 1. (Color online) Plot of the trace distance 
D(ps(t), Ps(t)) as a function of time, where p§(£) and p|(t) 
have been determined from Eq. (271. In both figures the 
horizontal line marks the upper bound of Eq. ( |17| ), and 
A = 0.1,g = 1 (arbitrary units), (a) Dynamics for two prod- 
uct total initial states which differ onl y by the environmental 
states and are given by Eq. (301 and ( |3l[ ) with |ai| 2 = 7/9, 



Q2| 2 = 8/9 and n — 7. (b) The two reduced states Ps{t) 
and p|(i) are obtained from the total initial states given 
by Eqs. ( 32 ) and ( 33 1 which have the same environmental 



marginal state, but different r educ ed syst em states and cor- 
relations. Parameters: a — i-^/3/7, /3 = \/4/7 and n = 1. 



while the right-hand side of Eq. ( 19 1 becomes 



D (p'seW, Ps(0)®Pe(0)) = \a(3\ 2 + \aP\. (35) 



Taking into account Eq. (23) and Eq. (24), for n,m ^> 
A 2 /4g 2 Eq. (34) explicitly reads 



D{p%{t),p 2 s (t)) = \a/3\ 2 [cos 2 (gy/m + It) - cos 2 (g^/nt) 
+ cos 2 (gy/rht) - cos 2 (gVn + It)] 
—S m , n -i Im {a*/3} sin (2gy/nt) , (36) 



which is an almost periodic function |39j since it repre- 
sents a linear combination of sine and cosine functions 
with incommensurable periods. The supremum of the 
attained values [ID] is less than or equal to 2|a/3| 2 if 
m 7^ n and m ^ n — 1, and equal to \a/3\ 2 + | Im{a*/3} | 
if m = n — 1. Thus, the inequality in Eq. ( |19[ ) is tight 
only for those initial states for which m = n — 1 and 
Re {a* j3} = (indeed, we have |a/3| 2 + \ap\ = 2|a/3| 2 if 

= 0). The special role of the 
1 can be traced back to the 



and only if either a = or , 
initial states with m = n - 
structure of the full unitary evolution given by Eq. ( 22 ) 



and to the presence of the creation and annihilation op- 
erators in the off-diagonal matrix elements. Their action 
generates, in fact, the last term in the modulus on the 
right-hand side of Eq. ( 36 ), which for m = n — 1 is neces- 
sary to reach the bound. If the relation n,m ^> A 2 /4<? 2 
is not satisfied, the supremum lies in general strictly be- 
low the bound even if the above mentioned conditions 
are fulfilled. This is a consequence of the fact that the 
periodic functions |c„(£)| 2 are then strictly less than 1. 
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III. GIBBS INITIAL STATE AND DYNAMICS 
OF THE TRACE DISTANCE 

We now extend our considerations to the evolution of 
the trace distance between a mixed correlated initial state 
and the tensor product of its marginals. Specifically, we 
will analyze the inequality given in Eq. ( 19 1 when the 



correlated initial state pse is the invariant Gibbs (ther- 
mal equilibrium) state corresponding to the full Hamil- 
tonian H of the model. For simplicity we will omit in 
the following the time argument zero. We first ana- 
lyze the total amount of correlations in the initial state 
D (pse, Ps ® Pe), i-e., the upper bound for the trace dis- 



tance according to Eq. (19). As we shall show below the 
main features of this bound can be explained in terms of 
the correlations in the ground state of the Hamiltonian 
H. We further study the behavior of the actual dynamics 
of the trace distance, which will turn out to reflect the 
characteristic features of the correlations in the Gibbs 
state. 



A. Correlations in the Gibbs state 

We consider the total initial Gibbs state 

1 



Pse 



(37) 



where H is the total Hamiltonian of the system given 
by Eq. (20), Z = Tre~^ H denotes the partition func- 
tion and /? = 1/kbT with kb the Boltzmann constant 
and T the temperature. To calculate the marginal states 
p s = Tr B e~ 0H jZ and p E = Tr s e~ pH /Z it is useful to 
obtain the matrix elements of pse with respect to the 
basis {\a, n)} already introduced in Sec. II B This can 



be done using the dressed states [41], i.e., the eigenvec- 
tors of the Hamiltonian H. These eigenvectors can be 
written as 



with 



|*+) =a„|l,n-l) + b n \0,n), 
l$n) = -b n \l,n- 1) +a n |0,ra) 
l$o> = |0,0>, 

1, 2, 3, . . . and 



(38) 



= \ I — , K — — , (39) 



where Q n = \/A 2 + 4g 2 n (see Eq. |24|). The corre- 
sponding eigenvalues are given by 



K„ = nuj + — ± 



2 ' 



Eo = 0. 



(40) 



Inverting Eqs. (38) with the help of the relations 

|0,n) =M$+)+a n |$->, 

|l,n) = a n+ i|$+ +1 > - b n+1 \$- +1 ), (41) 



one obtains the expressions 

.ran 

Poo ~z< 

pit = i< 



(e-^bl 



e ^ " a n I ■ 



pe: 



Pw l 



mn 

Poi 



+ 1 a. 



n+l 



S n +l,m ( e PE " +1 



e 

-PE. 



PE' 



a n+1 b n+1 , 
(42) 

which represent the matrix elements of the Gibbs state, 
Pse= f%p\<*,n){/3,m\- (43) 

a,/3,n,m 



Using this result together with Eq. ( 28 ) we see that the 
reduced system state is diagonal in the basis \a,n) and 
that the diagonal elements are given by p\\ = 1 — poo & n d 



Poo = 



n=0 



-PEi b 2 +£ -PE- 



(44) 



The reduced state of the environment is also diagonal 
since p nm = for n ^ m, and the diagonal elements can 
be expressed as 



p nn = 



l(e-^bl + e~^at 



Z 

+e 



-pe: 



+ 1 a. 



n+l 



(45) 



The product state constructed from the marginals is ac- 
cordingly of the form 



Ps® Pe = ^2p aa p nn \a,n)(a,n\. 



(46) 



Finally, the normalization constant Z can be written as 
Z = J2(e-^bl 



+ e -^+i<4 + 1 + e -^b 



2 

n+l 



(47) 



Starting from the above relations we can analyt- 
ically calculate the total amount of correlations of 
the Gibbs state, i.e., the quantity D{pse,Ps ® Pe)- 
To this end, we order the elements of the basis as 
{|0,0),|1,0),|0,1),|1,1),|0,2),|1,2),...}. The difference 
X = pse — Ps® Pe between the Gibbs state and its cor- 
responding product state can then be written in block 
diagonal form, 

/ D§ \ 

Dl p% 

pl\ Dl 

D\ pit 

pl\ Dl 



X 

















n,n+l 

Pi6 





n+l,n 

Poi 


Tjn+l 











'■■) 



(48) 
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where 



n,n 
ra.a 



(49) 



It is easy to demonstrate that D" = —Dfi, implying that 
the matrix of Eq. (48 1 has zero trace, as it should be. 



The eigenvalues of this matrix are simply given by the 
eigenvalues of the 2x2 block matrices plus the top left 
element D®. Hence, the total amount of correlations in 
the Gibbs state is given by 

D{pse,Ps®Pe) = \\Dl\ 



1 oo 

i: 

1 oo 

n=0 



n+l 



n,n+l 



D? + - d (D? - D^ 1 ) 2 + 4 for 1 



(50) 



This quantity depends on the model parameters Li, A 
and g which characterize the Hamiltonian described by 
Eq. (20), as well as on the temperature. In the fol- 



lowing we will focus in particular on the dependence of 
D{psEi Ps ® Pe) on the coupling constant g and on the 
inverse temperature /3 for fixed values of the other two pa- 
rameters (indeed from the expression of the Gibbs state 
it immediately appears that the dependence on one of 
the parameters can be reabsorbed into the others). 



B. Dependence on the ground state 



The behavior of the trace distance given by Eq. ( 50 ) as 
a function of (3 and g is plotted in Fig. |2]). We clearly see 
a non-monotonic behavior of the trace distance as a func- 
tion of both parameters. Focusing on the dependence on 
/3 for a fixed value of g, we observe that there is a sudden 
transition between two different kinds of behavior: Below 
a critical value of the coupling constant g, the trace dis- 
tance as a function of /3 exhibits an initial peak and then 
goes down to zero, see also Fig. (|3]a); above this criti- 
cal g it keeps growing to an asymptotic value different 
from zero, which we will discuss later on, as can be seen 
from Fig. ([3jb). On the other hand, the dependence of 
the trace distance on g for a fixed value of f3 shows some 
oscillations after a sudden growth which occurs at the 
critical g, see Figs. (2| and ([3]d). Quite remarkably, this 
means that the total amount of correlations of the Gibbs 
state can decrease with increasing coupling constant, as 
clearly observed in Fig. ([3jd). 

The above features can be explained considering that 
the trace distance D{pse, Ps ® Pe) quantifies the corre- 
lations of the Gibbs state pse and that the limit f3 — > oo 
corresponds to the limit of zero temperature, where the 
Gibbs state reduces to the ground state of the Hamilto- 
nian H. If all the eigenvalues given by Eq. (40) are non- 
negative the ground state is |$q ) = |0, 0) with eigenvalue 
zero. Of course, this is a product state and, therefore, the 




FIG. 2. (Color online) Plot of the correlations of the Gibbs 
state ( |37[ ) as a function of the inverse temperature /3 and of 
the coupling constant g according to Eq. (50 1. Parameters: 
ui — 3 and A = 0.5. 





FIG. 3. (Color online) (a, b, c, d From top left to bottom 
right) Sections of the plot in Fig. |2]) corresponding to g = 1.7, 
g = 5.5 , P = 5 and j3 = 8, respectively. The critical value of 
g is given by gi = 3.24, see Eq. (|51[). 



correlations of the Gibbs state approach zero for /? — > oo. 
This is what happens below the critical g. However, ac- 
cording to the level crossing described in Fig. Q, the 
Hamiltonian given by Eq. (20) has negative eigenvalues 



for larger values of the coupling constant g. In fact, it is 
easy to see from Eq. ( 40 1 that if 



9 > 9i 



-loA 



(51) 



then E^ < and, therefore, |0, 0) is no longer the ground 
state. Thus, we can then identify gi as the previously 
mentioned critical value of g, since for larger values the 
lowest energy state is I^T) which is an entangled state 
according to Eq. (38) with correlations a\b\ + a\bx dif- 
ferent from zero. But looking at the dependence of the 
different eigenvalues E~ on the coupling constant g, see 
Fig. Q, we can see that there is another critical point, 
let us call it g~2, where £^(32) = -ETfe) and after 
which E2 (g) < E^(g), i.e., \^) becomes the lowest 
energy state. We then have another value gs for which 
E3 (93) = -^(53)' so that for stronger couplings |$g) 
becomes the new ground state, and so on. Between two 
successive critical values g~i and the ground state of 
the Hamiltonian is |$~), whose correlations according to 
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Eq. ( 35 ) are given by 



D(p S E,Ps® Pe) = ajbj+ 



A 2 + V I/ A 2 



V 

This expression characterizes the asymptotic value of the 
correlations in the Gibbs state for /? — > oo and for g 
between g~i and g~i+i- We note that D(p$E,Ps <8> Pe) 
approaches the value | if we also let g — !• oo . As is shown 
in the Appendix, this asymptotic value corresponds in 
fact to the maximal possible value of the correlations for 
the present model. 



,e\ E 2 E iE - 



D(Pse.Ps®P£> 




FIG. 4. (Color online) (Top) Plot of the first energy eigen- 



values E^ , E2 , , E^ given by Eq. (401 as functions of g, 
Eg coincides with the :r-axis. (Bottom) Plot of the correla- 
tions of the Gibbs state as a function of g for j3 — 100, i.e., 
for approximately zero temperature; the other values are the 
same as in Fig. pi). The critical values of the correlations as a 
function of g exactly correspond to the level crossing points: 
When Eq = E^ there is a sudden increase and at the subse- 
quent points the dips occur. For this value of /3 the behavior 
described by the exact expression is well approximated by 



Eq. ( 52 1 between the dips and by Eq. ( 54 1 at the dips 



We see from Fig. Q that for small temperatures the 
correlations in the Gibbs state exhibit a dip at every gi 
with i > 1. Again, this feature can be explained con- 
sidering the ground level of the Hamiltonian given by 
Eq. (20 1. For g = §i the eigenspace of the lowest energy 
level is two-fold degenerate since E~(g~i) — E^^cn) and 



the Gibbs state reduces to 
1 



(i*r-i><*r-ii + i*r><*<"Q 



(53) 



where, again we have ordered the elements of the ba- 
sis as {|l,i- 2), |0,i- 1), 1), |0,z)}. Equation (ph 



can be directly obtained from Eq. (42), observing that 
for p — > 00 the only non- negligible terms are those in- 
volving the exponentials of (3E~_ 1 or (3E~ . Calculating 



now the corresponding product state and proceeding as 
done to obtain Eq. (50), or directly taking the limit of 
this equation for j3 — > 00 and g — g~i, one finds an explicit 
expression for the correlations of the mixed state given 
by Eq. p|: 



D {PsE,Ps® Pe) 



V 2 l7i + <*i + -M 2 +4e 2 | 
+~l7i +81- ]/ hi - Si) 2 + iej\ 



2 

+ 2 172 
1 



\fil2 



6 2 ) 2 



■4e 2 | 



+ d72 + h - ^(72-<5 2 ) 2 +4e 2 | + X 



(54) 



where 



b 2 b 2 



a = K 2 -i + a i) ; 71= 

*i = %■ - \ ("hi + «, 2 ) (ati + \ ^ = 
l2= b l-\(bU+bi)(^ 1+ b 2 )- £2 = -^; 

62 = f - I (aU + a 2 ) ; X = | {b 2 ^ + b 2 ) . (55) 



From the explicit evaluation of Eq. (|52j) and Eq. (54) 
for the different values of i, one can see that indeed the 
total amount of correlations of the mixed state given by 
Eq. ( 53 1 is smaller than the correlations of the dressed 
statesf^^) and |$~) giving its decomposition, which 
explains the emergence of the dips. Note however that 
the correlation measure given by D(pse, Ps ® Pe) is not 
a convex function on the space of physical states. 

The above arguments are summarized in Fig. Q. 
They explain the behavior of the correlations in the Gibbs 
state for small temperatures, i.e., for j3 — > 00. The effect 
of finite temperatures is to smoothen the dependence on 
0, as can be seen in Fig. @, gd) and |[3jc), such that 
the sudden increase at g — g~\ is less sharp and that the 
subsequent dips turn into oscillations which are more and 
more suppressed as the temperature increases. This be- 
havior is due to the fact that at finite temperature the 
Gibbs state has a non- vanishing admixture of |$7) for 
values of g which are smaller than g~i and, hence, the in- 
crease of the correlations starts before g = gi and is less 
sharp, as can be seen from Figs. Q and ([3]d). More- 
over, as a consequence of finite temperatures, the Gibbs 
state is a mixed state even between the critical values 
such that its correlations become smaller than in the 
zero temperature limit, which leads to a suppression of 
the oscillations. 



C. Time evolution of the trace distance 

The analysis performed so far concerns the correlations 
of the initial Gibbs state, i.e., the upper bound of the 
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trace distance between the reduced state Ps(t), evolving 
from an initial total Gibbs state, and the reduced state 
Pg(t), evolving from the corresponding product state, ac- 
cording to Eq. (19). We will now investigate the dy- 



namics of the trace distance D(pg(t), p 2 s {t)) and analyze, 
in particular, the dependence of the supremum of this 
function on the coupling constant and the temperature. 
As discussed before (see Sec. II A 2), the behavior of the 
trace distance between Pg(t) and pg(t) expresses the ef- 
fect of the initial correlations in the resulting dynamics. 
Moreover, its supremum as a function of time quantifies 
the amount of information which could not be initially 
retrieved by measurements on the reduced system only, 
but becomes accessible in the subsequent dynamics, thus 
making the two reduced states pg(t) and pg(t) more dis- 
tinguishable. 

Ta kin g as initial state p l SE the Gibbs state given by 
Eq. (37) and p 2 SE as the corresponding product of its 
marginals, we have p\{t) — Pg(0) since the Gibb state 
is invariant under the time evolution, and pg(0) = p|(0) 
because the corresponding open system initial states are 
identical. Thus, exploiting Eq. (45 1 we obtain the follow- 



ing explicit expression for the trace distance, 



DG£(*).p|(<)) 



(A 



oo 



+Poo ^2 np n 



(n + l)p n 



(56) 



For fixed values of the parameters characterizing the 
dynamics this expression describes a superposition of 
periodic functions with incommensurable periods, i.e., 
an almost periodic function as already encountered in 
Sec. |II B| An example for the trace distance dynamics is 
shown in Fig. The trace distance starts growing al- 
ready at the initial time and further oscillates with time, 
according to the almost periodic behavior described by 
Eq. J56l). 



As mentioned already the time dependence of the trace 
distance is solely due to the time evolution of the prod- 
uct state constructed from the marginals of the Gibbs 
state since the latter is invariant under the dynamics. 
It is the comparison between the two different reduced 
system states, namely between the states p$(t) = /9g(0) 
and Pg(t), which allows to obtain information initially 
not accessible with measurement on the reduced system 
only, and which enables the detection of correlations in 
the initial Gibbs state. 

The supremum of the trace distance in Fig. |5| is 
substantially smaller than the corresponding bound of 
Eq. (19). For large values of f3 and g the supremum 
can be estimated as follows. If the temperature goes 
to zero the Gibbs state approaches the projection onto 
the ground state which is given by I^X^ | for a fixed 
k, depending on the value of the coupling constant g. 
We suppose that g is different from the critical values 
cji. This implies p 00 = a 2 , pu — b 2 , together with 
p mm = Sfn^ai + Sm.k-ib 2 ,- For large values of g, which 



D(p>(t),pJ(t)) 

0.7 f 
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FIG. 5. (Color online) The trace distance D(p 1 s {t), p 2 s (t)) as 
a function of time according to Eq. ( |56| ; p's(t) is the state of 
the reduced system at time t obtained from an init ial total 
state Pse, where p\ E is the Gibbs state given by Eq. |37[ ) and 
p 2 SE is the corresponding product state. The upper horizon- 
tal line represents the bound given by the right-hand side of 
the inequality ( 19 \ which has been determined by Eq. ( 50 1 . 
Parameters: cj 



3, A = 0.5, g = 6 and /3 = 5. 



implies large values of k, we have ~ b^ ~ l/v2- Em- 
ploying further Eq. ( 23 ) , one thus obtains the estimate 



£>04(i),p|(i)) « i S m(2Vkgt)sm(gt/Vk) . (57) 

This shows that for large /3 and g the trace distance is 
bounded from above by \. 




FIG. 6. (Color online) The supremum of D(fi\(t), pg(t)) as 
a function of time versus the coupling constant g and the 
inverse temperature /3; pg(i) is obtained from an initial to- 
tal Gibbs state, p%{t) from the corresponding pro duct state, 
D{fis(f), P%(t)) is calculated according to Eq. (56 1. Parame- 
ters: ui — 3, A = 0.5. 

Figure ^ shows how the supremum of D{p 1 s {t) 1 p 2 g{t)) 
behaves as a function of the coupling constant g and the 
inverse temperature /3, keeping fixed oj and A. Exactly 
as for the correlations of the Gibbs state [compare with 
Fig. Q], we observe two qualitatively different kinds of 
behavior as a function of (3, for a fixed value of g. Be- 
low a critical g the supremum of the trace distance passes 
through maximum and then tends to zero; above the crit- 
ical value it tends monotonically to an asymptotic value 
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FIG. 7. (Color online) (a, b, c, d From top left to bottom 
right) The same as Fig. Q but for parameters g = 1.7, g = 
5.5, P = 5 and /3 = 100, respectively. For f} = 100, i.e., 
approximately zero temperature, the dips occur at the same 
values as the corresponding dips of the bound, see Fig. Q. 
For the case of finite temperature the dips are not suppressed, 
but they are shifted towards larger values of g. 



which is close to the estimate of \ determined above 



as 



illustrated in Figs. |7]a) and ([7jb). Moreover, considering 
the supremum of the trace distance as a function of g for 
fixed /3, after a sudden growth at the first critical g it ex- 
hibits some oscillations analogous to those of the bound. 
Comparing Figs. Q and we see that in the limit of 
zero temperature the bound and the true supremum of 
the trace distance both show a sudden increase and sub- 
sequent dips at the same values of the coupling constant 
g. This behavior can be explained by recalling the de- 
pendence of the energy spectrum of the Hamiltonian as a 
function of g in Fig. Q. At zero temperature the Gibbs 
state reduces to the ground level of the Hamiltonian. The 
discontinuous change in the ground level with varying g, 
i.e. the transition from to + |$ i+1 ))/\/2, im- 
plies a discontinuous change in the bound as well as in the 
supremum of the trace distance, thus leading to the dips 
appearing in Fig. j4| and Fig. Q. In fact, apart from 
fixing the bound at t le r.h.s. of Eq. ( [19] ), the Gibbs state 
determines both reduced states pg(t]and p s (t), arising 
from the initial total states pse and ps®PE respectively. 
Relying on Eq. ( 27 ) one can see that for A — the supre- 



mum of the trace distance is simply given by 1/4 for an 
initial correlated state pse = for any i > 0. 

This means that for zero detuning the supremum of the 
trace distance dynamics as a function of g at zero tem- 
perature takes the constant value 1/4, except at g = gi 
where the dips occur. The effect of a finite temperatures 
is slightly different for the bound and the supremum of 
the trace distance dynamics: With growing temperature 
the dips of the bound turn into oscillations which are 
more and more suppressed, but they occur at the same 
values of g. On the contrary, the dips of the true supre- 
mum, and its sudden increase as well, are not suppressed, 
but do change position, occurring at larger values of g. 



IV. CONCLUSIONS 



We have studied the influence of initial correlations be- 
tween system and reservoir on the dynamics of an open 
quantum system by means of the trace distance, consider- 
ing the paradigmatic and exactly solvable model provided 
by the Jaynes-Cummings Hamiltonian. First, we have 
analyzed the amount of correlations in the Gibbs state 
Pse as it is quantified by D(pse, Ps® Pe), where ps®pE 
denotes the product state arising from the marginals of 
Pse- The exact analytical expression of the latter quan- 
tity describes a non-monotonic behavior as a function of 
both the coupling constant and the temperature charac- 
terizing the dynamics. The same behavior is found for the 
supremum of the trace distance between the open system 
states p\{t) and p s (t) which evolve from pse and Ps®Pe, 
respectively. This enabled us to establish a clear con- 
nection between the correlation properties of the Gibbs 
state and basic features of the subsequent open system 
dynamics, in particular, the amount of information which 
is initially inaccessible for the open system and which is 
uncovered during its time evolution. The dynamical be- 
havior of the trace distance between the reduced system 
states p s (t) and p s (t) thus provides a witness for the 
correlations of the total system's initial state. 

As we have shown for the case at hand, at zero temper- 
ature sudden changes in the supremum over time of the 
trace distance can be traced back uniquely to discontinu- 
ous changes in the structure of the total system's ground 
state and to its degree of entanglement which, in turn, is 
caused by crossings of the energy levels of the total sys- 
tem Hamiltonian. It is important to remark that, as we 
have demonstrated, clear signatures of these discontinu- 
ities are still present at finite temperatures. Note that to 
reconstruct the trace distance dynamics, in order to de- 
tect correlation properties of the ground state, one only 
needs to follow the evolution of the open system state 
obtained from the initial total product state ps <8> Pe- 



The bound given by the right-hand side of Eq. ( 19 ) is 



able to represent qualitatively the non-trivial behavior of 
the maximum of the trace distance between p s (t) and 
Ps(t), as a function of the different parameters charac- 
terizing the Hamiltonian and the temperature. While for 
the sudden transition between the two different asymp- 
totic regimes as a function of j3 it is clear that the effective 
maximum of D(p 1 s (t), p s (t)) has to reproduce the behav- 
ior of the bound, it is quite remarkable that also in the 
second case, where the bound is sensibly different from 
the effective maximum and from zero, both these quanti- 
ties show an analogous non-monotonic behavior. Finally, 
we note that it must be expected that the general features 
found here for the correlated Gibbs state hold true also 
for other correlated initial states, e.g., for correlated non- 
equilibrium stationary states, as long as the latter involve 
discontinuous, qualitative changes under the variation of 
some system parameters. 
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Appendix: General bound for the correlations of a 
quantum state 

Throughout the article we have used the quantity 
D{psEi Ps®Pe) as a measure for the total amount of cor- 
relations of given state pse- On the ground of extensive 
numerical simulations we conjecture that this quantity 
satisfies the inequality 



D(pse,Ps <8> Pe) < 1 



1 

N 2 ' 



(A.l) 



where N denotes the minimum of the dimensions of %g 
and %e- For the example studied in this paper we have 
N = 2 and, hence, D(p SE ,ps ® Pe) < f • 

To our knowledge there exists no general mathematical 
proof for the inequality (A.l). However, one can easily 



prove that this inequality is saturated if pse — IV , )('0l 
is a pure, maximally entangled state. To show this we 
first note that for a maximally entangled state vector \ip) 
the marginal states are given by ps — Ps/N and pe — 
Pe/N, where Ps and Pe are the projections onto the 
subspaces of fis and He, respectively, which are spanned 
by the local Schmidt basis vectors with nonzero Schmidt 
coefficients. Hence, D(p$E, Ps ® Pe) is given by \ times 
the sum of the absolute eigenvalues of the operator 



X = \ip)(ip\ 



1 

N 2 



Ps ® Pe- 



(A.2) 



Obviously, \ip) is an eigenvector of X corresponding to 
the eigenvalue 1 — 1/N 2 . Moreover, all vectors which 
are perpendicular to \ip) and belong to the support of 
Ps®Pe are eigenvectors of X with the eigenvalue —1/N 2 . 
Thus, X has one non-degenerate eigenvalue 1 — 1/N 2 , and 
one eigenvalue —1/N 2 which is (N 2 — l)-fold degenerate, 
while all other eigenvalues of X are zero. Therefore we 
have 



D(pse,Ps®Pe) = 



1 



1- ^ + (^ 2 -l)^ 
N 2 N 2 



= 1- 



1 



which proves the claim. 



N 2 
(A.3) 
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